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The adsorption of two different organic molecules cyanidin glucoside (C 21 O 11 H 20 ) and TA-St-CA 
on anatase (101) and (001) nanowires have been investigated using the standard and the range sep¬ 
arated hybrid density functional theory calculations. The electronic structures and optical spectra 
of resulting dye-nanowire combined systems show distinct features for these types of photochro- 
mophores. The lowest unoccupied molecular orbital of the natural dye cyanidin glucoside is located 
below the conduction band of the semiconductor while, in the case of TA-St-CA, it resonates with 
the states inside the conduction band. The wide-bandgap anatase nanowires can be functionalized 
for solar cells through electron-hole generation and subsequent charge injection by these dye sen¬ 
sitizers. The intermolecular charge transfer character of Donor-7r-Acceptor type dye TA-St-CA is 
substantially modified by its adsorption on Ti02 surfaces. Cyanidin glucoside exhibits relatively 
stronger anchoring on the nanowires through its hydroxyl groups. The atomic structures of dye- 
nanowire systems re-optimized with the inclusion of nonlinear solvation effects showed that the 
binding strengths of both dyes remain moderate even in ionic solutions. 

PACS numbers: 71.15.Mb, 68.47.Gh 


I. INTRODUCTION 

The wide bandgap metal oxide, Ti02 has gained an 
increased attention since the discovery of its ability to 
carry out hydrolysis under UV irradiation by Fujishima 
and Honda. In addition to showing such an excellent 
photocatalytic performance, Ti 02 has also become the 
material of choice as the anode in dye sensitized solar cell 
(DSSC) applications due to its favorable electrochemical 
and charge carrier conduction properties. [HHin] 

In terms of more efficient utilization, titania pos¬ 
sesses particularly importance with its nanocrystalline 
forms. [nHn Quasi-one-dimensional Ti02 nanostruc¬ 
tures offer high surface-to-volume ratios that is desirable 
to improve efficiencies of photovoltaic and photocatalytic 
processes. [TMT5] 

The (001) and (101) terminations of the anatase poly¬ 
morph are known to exhibit remarkably higher photo¬ 
catalytic activity relative to the surfaces of the rutile 
phase. [inMl] Moreover, the naturally occurring anatase 
form has been reported to be the most stable phase of 
Ti 02 in nanodimensions. [HHMj 

In a basic DSSG operation, many processes take place 
influencing overall device performances. These are the 
electron-hole generation upon visible light absorption by 
the dye sensitizer, the charge carrier injection from the 
dye to the conduction band (GB) of the oxide electrode, 
regeneration of the ground state of the dye by a redox re¬ 
action through liquid iodide/triiodide electrolytes. [51 ITU] 
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Here one of the main concerns is the reduction of pho¬ 
togenerated electron-hole recombination rate. Another 
point is the binding strength of the sensitizer molecule 
to the oxide and the reliability of this system in the ionic 
solvent. Therefore, the type of the dye molecule and the 
surface properties of the semiconductor plays a key role 
in the optimization of such processes. 

Recent experiments have shown that one-dimensional 
nanostructures have several advantages over nanopartic¬ 
ulate Ti02 films. [T21 [TK] First of all they have higher 
surface-to-volume ratios allowing many active sites to 
come into contact with light harvesting molecules. Sec¬ 
ondly, their one-dimensionality forms a natural pathway 
for the charge carrier conduction from the point of in¬ 
jection to the anode. Moreover, nanowires or nanotubes 
can exhibit a band transport rather than a hoping mech¬ 
anism between nanoparticles. Another crucial aspect is 
that the electron-hole recombination rates were observed 
to be ten times lower in nanotube-based films in compar¬ 
ison with films made from nanoparticles. m 

One of the key issues in DSSG is the type of the sen¬ 
sitizer dye and its interaction with the Ti 02 nanostruc¬ 
tures. Commercially available Ru-based molecular com¬ 
plexes have shown up to 11% solar energy conversion 
efficiencies. [HIIMHSO] These metal-driven dyes can per¬ 
form spatial charge separation and fast injection rates. 
Meanwhile, researchers have also focused on finding nat¬ 
ural alternatives. Recently, antenna type novel organic 
donor-TT-acceptor (D-tt-A) dyes have been proposed as 
sensitizers to achieve charge transfer excitations reduc¬ 
ing recombination rates that is the basic concern in sim¬ 
ple skeleton light harvesting molecules. For instance, tt- 
conjugated organic complexes with tetrahydroquinoline 
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FIG. 1. Optimized adsorption geometries of cyanidin glucoside on anatase (OOl)-nanowire and (lOl)-nanowire 


moiety as donor and a cyanoacrylic acid moiety as the 
acceptor were experimentally synthesized and theoreti¬ 
cally studied. pTTHHi] From this class, a highly efficient 
organic dye, TA-St-CA, contains a 7r-conjugated oligo- 
phenylenevinylene electron donor-acceptor moiety and 
a carboxyl group as anchoring group. |551 - HT| Hwang et 
al. achieved % 9.1 photo-to-electric conversion efficiency 
with an open circuit voltage of 743 mV by designing a 
low cost TA-St-CA based DSSC.[55] 

Natural dye pigments with a simple carbon skeleton 
structures such as cyanidin dyes are eco-friendly, widely 
available and cheap to produce. For instance, antho- 
cyanin (cyanidin-3-O-glucoside) can easily be extracted 
from plants. [42] The cyanidin family is well known and 
is proposed as an alternative to other dye sensitizers. m- 
[5^ However, reports indicate the overall solar-to-electric 
energy conversion efficiencies below the current require¬ 
ments. A deeper understanding of the bottlenecks in 
their DSSC applications is still needed. 

In this study, we theoretically investigated the adsorp¬ 
tion modes, electronic structures and absorption spectra 
of two different types of dyes on anatase Ti 02 nanowires 
having (001) and (101) facets. For this reason, we con¬ 
sidered the D-tt-A type organic complex TA-St-CA and 
the natural chromophore anthocyanin as dye sensitizers. 
We used both the standard and the screened exchange 


hybrid density functional theory (DFT) calculations to 
shed light on the main differences between the resulting 
nw-|-dye combined systems in vacuum as well as in solu¬ 
tion treated with a modern nonlinear polarizable contin¬ 
uum model (PCM). 


II. COMPUTATIONAL DETAILS 

The DFT calculations have been performed using the 
projector-augmented wave (PAW) methodjTO] as im¬ 
plemented in the Vienna ah-initio simulation package 
(vASP).|7ii izg Single particle electronic states have 
been expanded using plane wave basis sets up to a ki¬ 
netic energy cutoff value of 400 eV. We used the standard 
generalized gradient approximation (GGA) to describe 
the exchange-correlation (XG) effects with the semi¬ 
colon Perdew-Burke-Ernzerhof (PBE)|73] functional as 
well as the contemporary hybrid XG functional proposed 
by Heyd-Scuseria-Ernzerhof (HSEl. [7iH7B] The latter is 
a nonlocal, screened Coulomb potential scheme with a 
range-separation. HSE functional largely heals the in¬ 
herent bandgap underestimation of the standard DFT. 
To do so, hse[71H7S| proposed to partially admix the 
exact Fock and the PBE exchange energies in the short 
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FIG. 2. Optimized adsorption geometries of TA-St-CA on anatase (OOl)-nanowire and (lOl)-nanowire 


range (SR) part as 

ijHSE ^ ^ (1 _ „)i5PBE.SR(^) ^ 

where a is the mixing coefficient |79j and w is the range 
separation parameter. [ziHini Meanwhile, the long range 
(LR) part of the PBE exchange and the full PBE cor¬ 
relation energies are included in the HSE functional. In 
this way, the lack of proper self-interaction cancellation 
between the Hartree and exchange terms of the standard 
DET is partly avoided leading to a substantial correc¬ 
tion of the bandgap underestimation. In addition, the 
tendency of the standard DET to give overly delocal¬ 


ized charge density distributions is also corrected to some 
extent. Therefore, this range separated hybrid density 
functional approach not only improves the bandgap re¬ 
lated properties over the standard exchange-correlation 
(XC) schemes but also offers a better description of lo¬ 
calized states such as the Ti 3d states of Ti02 or iso¬ 
lated gap states caused by various impurity atoms. [II1[T3 
Hence, the HSE functional offers an improvement over 
PBE-calculated electronic and optical properties while 
the difference in the optimization of lattice structures is 
relatively less recognizable between PBE and HSE cal¬ 
culations. In fact, for dye-|-nanowire systems considered 














4 


in this work, both XC functionals gave similar relaxed 
geometries and binding modes. 

In order to obtain absorption spectra as the imagi¬ 
nary part of the dielectric function, £ 2 (w), from a density 
functional calculation, one considers the transitions from 
occupied to unoccupied states within the first Brillouin 
zone as the sum, 

4^/3 1™, ^ H 2wk^(eck - e«k - w) 

^ ('^ck+eQ.g ('^ck+e^g 

where the c and v show empty and filled states respec¬ 
tively, Uck are the cell periodic part of the orbitals and 
Wk are the weight factors at each k-r)oint.[HD] 

The nw(OOl) and nw(lOl) nanowire models are con¬ 
structed from the anatase Ti02 bulk structures. They 
are considered in large tetragonal supercells with dye ad¬ 
sorbates such that the periodicity of the cell is chosen to 
be five times larger than one unit cell length along the 
nanowire axis to isolate the molecules from each other 
on the nanowires (as in Fig. and Fig. |^. In order to 
avoid unphysical interactions between the periodic im¬ 
ages of the dye-fnw structures, the supercells contain at 
least 20 A of vacuum separations along both of the lattice 
translation vectors perpendicular to the one along the 
nanowire axis. We fully optimized initial geometries by 
minimizing the Hellman-Feynmann forces on each ionic 
core to be less than 0.01 eV/A based on the conjugate- 
gradients algorithm. None of the atoms were frozen to 
their bulk positions during these relaxation procedures. 
In this way, the bare nanowire models were previously 
shown to maintain the anatase structure without a ma¬ 
jor lattice distortion. [ST] 

The effect of the solvent environment on the elec¬ 
tronic structure of the dye-nanowire composed systems 
has been studied using the nonlinear polarizable contin¬ 
uum model (PCM) as implemented in the open-source 
code JDFTx.jSSHM] Since the nonlinear models incor¬ 
porate the dielectric saturation effect, the nonlinear [82] 
PCMs are expected to be more accurate than similar 
linear models in the presence of highly polar structures 
such as Ti02. In such models, the contribution to the 
dielectric function arising from the rotations of solvent 
molecules are modeled as a field of interacting dipoles, 
whose response function saturates with increasing ex¬ 
ternal field. The cavity surrounding the solute is con¬ 
structed self-consistently from its electron density, where 
the dielectric function of the solvent turns on at a criti¬ 
cal electron density contour. The numerical value of the 
critical electron density and the effective surface tension 
of the solute-solvent interface are highly solvent depen¬ 
dent. For solvents composed of small and highly polar 
molecules (such as water), the effective surface tension at 
the interface is positive; whereas for solvents with large 
molecules and strong dispersion (van der Waals) inter¬ 
actions (such as chloroform) this effective tension often 
has a negative sign. Additional technical details on the 


polarizable continuum model used as well as the numer¬ 
ical values for the critical electron densities and effective 
surface tensions of the solute-solvent interfaces can be 
found in the relevant publications. [5^ IBTl 155] 


III. RESULTS & DISCUSSION 

The natural dye (cyanidin 3-0-glucoside) has been con¬ 
sidered as a sensitizer for DSSC applications by various 
groups. [HHSni The main advantages over the other or¬ 
ganic dyes have been pointed out as the simple and low 
cost production together with its wide availability in na¬ 
ture. Theory based studies on this class of dyes are rather 
rare and mostly limited to isolated molecules in the gas 
phase. We considered various probable initial adsorp¬ 
tion configurations of the cyanidin dye on both (001) and 
(101) anatase nanowire models, we carried out the full 
optimization of the atomic coordinates using both PBE 
and HSE functionals. Results indicate energetically fa¬ 
vorable two different binding modes on both nanowires as 
shown in Fig. [^ In order to differentiate between these 
two modes we refer them as the physical and chemical 
bindings. In the physical binding, the tail oxygen of the 
cyanidin part and the oxygen of the OH group at the 
end of glucoside moiety interacts with two five-fold coor¬ 
dinated surface Ti atoms (see Fig[^ and Fig[^). In the 
chemical binding mode, the OH group additionally loses 
its hydrogen to the nearest surface oxygen site on both 
nw(OOl) and nw(lOl) models as shown in Fig and 
Fig[T}i, respectively. Hence, the calculated adsorption 
energies presented in Table |l] indicate stronger binding 
for the latter one. 

The HSE-calculated Ti-0 bond lengths between the 
cyanidin dye and nw(OOl) are 2.04 A (on the cyanidin 
part) and 2.22 A (on the glucoside side) for the physical 
adsorption mode. The second bond shortens to 1.88 A in 
the chemical binding case. On nw(lOl), the correspond¬ 
ing bonds are 2.35 A and 1.99 A for the physical binding 
while they are 1.97 A and 1.91 A for the chemical binding 
mode. 

The geometry optimizations using the standard PBE 
and modern HSE functionals gave similar final structures 
except the chemical binding mode of cyanidin molecule 
on the (101) nanowire. PBE calculations predict much 
stronger interaction between the natural dye and the 
(101) surface of the oxide. The PBE functional leads to 
considerable local distortion around the adsorption site 
where the bond lengths abruptly change. For instance, 
PBE-relaxation breaks the bond between the surface Ti 
which the dye is anchored at and the adjacent surface 
O which captures the hydrogen from the dye. The sep¬ 
aration between them increases from the typical surface 
0-Ti bond length of 1.97 A to 3.81 A. While PBE causes 
such a local reconstruction, HSE calculations yield an 
adsorption geometry where the (101) nanowire keeps its 
bare surface structure. Therefore, calculated binding en¬ 
ergies for this case are 2.35 eV and 1.20 eV with PBE 
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FIG. 3. Charge densities of the highest occupied states (on 
the left) and the lowest unoccupied states (on the right) of 
cyanidin glucoside dye (top row), dye+nw(101) (middle row), 
and dye+nw(001) (bottom row). 


and HSE schemes, respectively (see Table]!]). Experimen¬ 
tal studies does not report a binding at the chemisorp¬ 
tion level to support the PBE predictions for this specific 
case. Therefore, this can be seen as one of the examples 
where the standard exchange-correlation functionals end 
up with peculiar results. 

The TA-St-CA molecule on the (001) and (101) facets 
of the anatase nanowires has two different low energy 
adsorption structures, the monodentate and the biden- 
tate binding modes as shown in Fig]^ The monodentate 
binding portrays perpendicular alignment with respect to 
the nanowire axis as a result of the single bond formation 
between the tail oxygen and surface Ti atom. The biden- 
tate mode is similar to the cyandin case because of the 
additional H transfer from the OH group at the tail to 
the nearest surface oxygen site on both of the nanowire 
types. The loss of hydrogen from the dye to the sur¬ 
face enables another 0-Ti bond formation between the 
molecule and the oxide. 

In the bidentate case, both 0-Ti bonds are 2.01 A be- 



FIG. 4. Gharge densities of the highest occupied states (on 
the left) and the lowest unoccupied states (on the right) of 
TA-St-GA dye (top row), dye-l-nw(lOl) (middle row), and 
dye-l-nw(OOl) (bottom row). 


tween TA-St-CA and nw(lOl). Likewise, they are 2.03 
A on nw(OOl). Monodentate TA-St-CA forms a single 
bond of 2.14 A on nw(lOl) while it is 2.36 A on nw(OOl). 
These results expectedly indicate stronger binding in fa¬ 
vor of the bidentate adsorption. A similar conclusion can 
be drawn from the calculated binding energies in Table ]^ 
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TABLE I. Calculated adsorption energies, -Bads, of the dye-nanowire systems 


Dye 



@( 001 ) 




@( 101 ) 


PBE 

HSE 

PBE-tPCM^ 

PBE+PCM^ 

PBE 

HSE 

PBE+PCM^ 

PBE+PCM^ 

cyanidin glucoside (physical) 

-1.35 

- 1.22 

-1.13 

-0.78 

-1.47 

-1.42 

-1.16 

-0.74 

cyanidin glucoside (chemical) 

-1.67 

-1.30 

-1.33 

-0.89 

-2.59 

- 1.20 

-2.08 

-1.44 

TA-St-CA (monodentate) 

-0.24 

-0.28 

-0.03 

0.18 

-0.74 

-0.72 

-0.50 

-0.26 

TA-St-CA (bidentate) 

-1.05 

- 1.02 

-0.81 

-0.58 

-0.84 

-0.93 

-0.64 

-0.43 


^ in CH 3 CI using nonlinear PCM 
^ in H 2 O using nonlinear PCM 


Molecular complexes with anchoring groups are rather 
obvious to yield larger adsorption. However, the main 
difference between TA-St-CA and cyanidin dyes is not 
their sizes. The most important factor is that while TA- 
St-CA is specifically designed to do an intramolecular 
charge separation from the antenna part to the accep¬ 
tor moiety which attaches to the oxide surface, both the 
HOMO and LUMO charge densities are localized on the 
cyanidin part of the cyanidin-3-O-gloucoside. The gluco- 
side group does not give any contribution to the frontier 
molecular orbitals in the case of cyanidin dye. This can 
be seen from the calculated charge density distributions 
for the corresponding states of the molecules in their gas 
phase as shown at the first rows of Fig and Fig. 

In fact, the comparison of the charge density redistri¬ 
bution features for the lowest lying optical excitation be¬ 
comes more important when these molecules are attached 
to the oxide surfaces. Therefore, the spatial charge den¬ 
sities of the highest occupied and the lowest unoccupied 
states of dye-l-nw combined systems have been calculated 
to discuss the charge injection features of these two dif¬ 
ferent types of dyes. The HSE results are presented in 
the middle and at the bottom rows of Fig. and Fig. 
for the nw(lOl) and the nw(OOl) cases, respectively. The 
charge densities of the frontier molecular orbitals of the 
cyanidin dye remain very similar to their gas phase dis¬ 
tributions even if it forms two 0-Ti bonds with both the 
anatase (001) and (101) surfaces. These two states also 
appear in the band gap of nanowires as well-localized 
isolated states, one being filled and the other one being 
empty (see Fig. [^. Hence, the lowest vertical excita¬ 
tion does not involve a charge injection to the CB of the 
semiconductor. Such an excitation is prone to yield an 
electron-hole recombination. Therefore, HSE-calculated 
charge density results might partially explain why this 
type of natural dye pigments end up with relatively low 
incident photon to current efficiencies (IPCE). 

The charge density distributions of the frontier orbitals 
of an isolated TA-St-CA molecule in vacuum have been 
investigated by Zhang et al. m and Mohammadi et al. [41] 
using hybrid DFT calculations. Their gas phase re¬ 
sults agree with our calculations for the isolated molecule 
case as shown in the top row of Fig. In the case of 


TA-St-CA on Ti02, the intramolecular charge transfer 
character from the triphenylamine part as donor to the 
cyanoacrylic acid group as the acceptor seems to be sig¬ 
nificantly altered after the molecule is adsorbed on the 
anatase nanowires. The lowest lying excitation involve a 
transition from the HOMO-like state which is spatially 
well-localized on the antenna moiety to the CB of the 
anatase as shown for the TA-St-CA-|-nw combined sys¬ 
tems in Fig. [^ Therefore, our calculations indicate that 
upon adsorption the loss of donor-7r-acceptor gas phase 
feature of TA-St-CA changes in favor of better charge 
injection into the CB of the oxide. This might also be 
interpreted as accounting for the reduction of recombi¬ 
nation rates in the case of phenylenevinylene-conjugated 
D-tt-A type sensitizers. 

The adsorption energies of the molecules have been 
calculated using 

Eads — Edye+ nw (Edye + ^nw) 

where Edye+nwj Edye, and are the relaxed super¬ 
cell energies of dye adsorbed nanowire, the isolated dye 
molecule, and the bare anatase nanowire, respectively. 
Both the PBE and HSE functionals were used for the 
vacuum calculations as presented in Table |T| Then, we 
obtained adsorption energies in solution for chloroform 
and water using a new non-linear PCM|52| with the PBE 
functional. Our tests show that similar conclusions can 
be drawn when the solvent effects are included within 
the hybrid HSE scheme. The vacuum results with PBE 
and HSE functionals are only slightly different from each 
other except the chemical binding of the cyanidin dye on 
the nw(lOl). For this case, the energy difference of 1.39 
eV is due to the fact that the PBE functional overesti¬ 
mates the bonding between the hydrogen atom captured 
from the dye with the nearest surface oxygen, which re¬ 
sults in an additional local modification of atomic posi¬ 
tions relative to the HSE-optimized structure where sur¬ 
face Ti-0 row seems to be disturbed less. 

As expected, single bond formation leads to a weaker 
adsorption energy for each type of dye. The loss of H 
from the OH group of the dye to the nearest surface oxy¬ 
gen causes the formation of a second bond. This sit¬ 
uation enhances the binding appreciably, especially in 
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FIG. 5. Densities of states for the dye+nanowire systems (in arbitrary units) calculated using the PBE and HSE functionals. 
The partial DOS contributions of the dye molecules are shown as dark (blue) shades. 


the bidentate mode of TA-St-CA. Clearly, the binding of 
the molecules appears to be noticeably stronger in vac¬ 
uum. As a nonpolar solvent with a dielectric constant of 
4.8, chloroform has a little effect on the adsorption ener¬ 
gies. However, water, being a polar solvent, weakens the 
bond(s) between the dye and the nanowire, considerably. 
For instance, the results show that monodentate TA-St- 
CA is washed away from the surface. In the other cases, 
the dye molecules can keep moderate binding with the 
oxide even in water environment. 

We note that the binding energies change in the pos¬ 
itive direction as seen in Table HI The reason is that the 
binding sites on the Ti 02 nanowire and the dye molecule 
interact strongly with the solvent. On the other hand, 
such an interaction is absent in a vacuum calculation 
leading to a more negative binding energy. When we 
check the differences between the calculated Kohn-Sham 
eigenvalues in vacuum and in solution, we see that the en¬ 
ergy levels are slightly shifted with respect to each other. 


Since H2O is a more polar solvent than CHCI3 with a 
higher dielectric constant, it interacts more strongly with 
the binding sites; we therefore see that binding energies 
in H 2 O are more positive than the binding energies in 
CHCI3. 

Since the dye molecules develop a bonding interaction 
with the oxide surface, dye-related HOMO-like levels ap¬ 
pear in the band gap of the oxide as isolated and well- 
localized states. In addition, their LUMO levels show a 
strong dispersion inside the CB of Ti02. Those can no 
longer be considered as molecular energy levels. There¬ 
fore an optical excitation starts from the HOMO-like dye- 
related state to the states in the CB. 

For an efficient DSSC, the lowest lying absorption 
peaks (in the visible range) mainly involve transitions 
from the HOMO-like dye-related states which appear 
above the VB to the states in the CB of Ti02. If the 
dye molecules had broken apart and been dissolved in 
solution one could have expected the loss of main fea- 































































































tures of the absorption peaks. However, we have shown 
that the water or chloroform can only weaken the bonds 
between the molecules and Ti02. When we check the 
differences between the calculated Kohn-Sham eigenval¬ 
ues in vacuum and in solution, we see that the single 
particle energy levels are slightly shifted with respect to 
each other. Therefore, from a theoretical point of view, 
the difference in the calculated electronic properties be¬ 
tween the standard PBE and modern HSE exchange- 
correlation schemes is much larger than the energy shifts 
due to PCM. 
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FIG. 6. First two rows display the calculated absorption spec¬ 
tra of the dye-|-nanowire combined systems : (a) physical (b) 
chemical binding of cyanidin glucoside on the anatase nw(OOl) 
and nw(lOl), respectively. The last two rows show the absorp¬ 
tion spectra for (c) monodentate TA-St-CA and (d) bidentate 
TA-St-CA on anatase nw(OOl) and nw(lOl), respectively. 


In order to discuss the electronic structures of the 
dye-|-nanowire combined systems the partial density of 
states has been calculated for each binding mode using 
both the PBE and HSE functionals as presented in Fig.[^ 
Since, the DOS structures calculated using the PBE and 
HSE functionals were aligned with respect to their deep 
core states, the HSE VB edge lies lower than the PBE 
one. Therefore, in the case of HSE, several new dye- 
related states fall in the band gap, which seem to occupy 
the top of the VB in the PBE DOS structures. The inclu¬ 
sion of partial short ranged exact exchange corrects the 
description of Ti Sd states reasonably. Since the CB of 
Ti02 is composed of these d bands, the CB gets shifted 
up to higher energies. Therefore, HSE functional heals 
the band gap underestimation of the standard exchange- 
correlation schemes not only by shifting up of the CB 
edge but also by lowering the VB edge. Meanwhile, the 
presence of new dye related gap states cause a signifi¬ 


cant narrowing of the band gap, which functionalizes the 
combined system to be active in the visible region. 

The absorption spectra of the dye-fnanowire systems 
have been obtained at the PBE and HSE levels by cal¬ 
culating the dipole transition matrix elements between 
the occupied and the empty states. Similar computa¬ 
tions were performed previously for the bare anatase 
nanowires. [M] The presence of new isolated and occu¬ 
pied gap states originating from the dye adsorbates just 
above the VB edge of the oxide is desirable for func¬ 
tionalization of Ti 02 in the visible part of the spectrum. 
When one compares the calculated optical spectra, simi¬ 
larities can be found between the two nanowire types as 
well as between the results of PBE and HSE function¬ 
als as shown in Fig. Among these characteristics, the 
first absorption peaks in each case are associated with 
the vertical transitions from the dye-related HOMO-like 
gap states to the lowest lying unoccupied states. Those 
sharp peaks indicate a significant red shift of the absorp¬ 
tion threshold into the visible region which is favorable 
for titania based photovoltaic applications. Although the 
overall features look alike, PBE-calculated spectra is also 
considerably red shifted relative to those of the HSE due 
to the local density approximation (EDA) to exchange- 
correlation effects giving rise to an underestimation of 
the band gap of Ti02. Moreover, the red shifting of the 
PBE spectra with respect to the HSE results gets even 
larger for the higher lying transitions. 

Experimentally, Senthil et a/. [42] reported the first 
absorption peak of cyanidin-3-O-glucoside adsorbed on 
Ti02 nanoparticles at 2.22 eV. Although the efficiency 
they have obtained with this dye is still below the require¬ 
ments as they said, the main reason remained unclear. 
In our calculation for the cyanidin glucoside on both of 
the nanowire types, the first absorption peak positions 
are found around 1.5 eV and 2.2 eV with PBE and HSE 
functionals, respectively, which coincide with the energy 
difference between the isolated dye-related gap states as 
seen in the top two panels of Fig. Therefore, the low¬ 
est lying transitions in each case essentially involve an 
intramolecular excitation from the HOMO-like state to 
the LUMO-like state in the band gap as seen in first two 
panels of Fig. The corresponding charge density dis¬ 
tributions can also be seen in Fig. Although the first 
absorption peak seems to drive the photoresponse of the 
combined system to visible region, the associated exci¬ 
tation takes place on the dye itself without any sign of 
charge injection into the semiconductor. This might be 
seen as one of the factors why these type of natural dyes 
have relatively low incident photon to current efficiency 
(IPCE) observed in experiments. |42j 

Recently, Hwang et al. reported an overall solar-to- 
energy conversion efficiency of 9.1 % with TA-St-CA 
sensitizer. |35j The experimental optical absorption peaks 
were found around 2.4 eV.|35l|36] This value agrees with 
our HSE results for the bidentate TA-St-CA. The re¬ 
markable difference in the efficiencies between the nat¬ 
ural cyanidin and organic TA-St-CA sensitizers needs an 
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understanding from a theoretical perspective. First of 
all, the adsorption of TA-St-CA causes the LUMO to 
resonate with the Ti 3d states and get delocalized in the 
CB of the anatase nanowires. Therefore, the first peaks 
corresponding to absorption in the visible region are due 
to transition from the dye-related HOMO-like gap states 
to the states in the CB as seen in the bottom two panels 
of Fig. The inclusion of screened exact exchange in 
the HSE functional shifts the CB edge to higher ener¬ 
gies such that the gap values become ~ 1.3 eV whereas 
the PBE gaps are considerably underestimated as seen 
in Eig. Consequently, the hrst peak positions agree 
with the energy difference of those states involved in the 
transition. When we look at the corresponding charge 
density plots in Fig. we see a charge injection into 
the oxide through electron-hole generation upon visible 
light absorption. In comparison to the cyanidin case, the 
appearance of a group of dye-related states just above 
the VB edge in the HSE results translates into an ab¬ 
sorption ability of the TA-St-CA-|-nanowire system in a 
wider range of the visible spectrum with varying oscilla¬ 
tor strengths. For the TA-St-CA sensitizer, transitions 
might arise from these group of gap states below the 
HOMO-like level to the states in the CB causing shal¬ 
low peaks around 3 eV shown in Fig. and d. The 
results indicate that the cyanidin dye is not successful 
in sensitizing anatase nanowires in the region between 
the first peaks and the main body of the higher en¬ 
ergy contributions. This is much more pronounced in 
the spectra obtained with HSE functional. Higher fre¬ 
quency contributions to absorption in the UV part of the 
spectrum are mostly associated with the interband tran¬ 
sitions from the states in the VB to the states in the CB 
of the dye-l-nanowire system. Both molecules bring the 
lowest lying peak which drops the absorption threshold 
into the visible region as seen in Fig. for both nanowire 
types. The TA-St-CA-|-nanowire combined system has 
more favorable optical properties than the cyanidin glu- 
coside-l-nanowire structure. 


IV. CONCLUSIONS 

The natural chromophore cyanidin-3-O-glucoside and 
the D-tt-A type organic dye TA-St-CA have been sep¬ 
arately considered on anatase Ti 02 nanowires having 
either (101) or (001) facets. The standard and the 
screened Coulomb hybrid density functional theory cal¬ 
culations were performed to understand the adsorption 
modes, electronic structures and absorption spectra of 
dye+nanowire combined systems. The binding energies 
have been obtained in vacuum and in solution using a 
new non-linear PCM. Our results indicate a significant 
band gap narrowing upon adsorption on both of the 
nanowire facets due to the appearance of a number of 
isolated dye-related gap states. Their number increases 
when a second bond is formed between the dye molecule 


and the oxide. The HSE-calculated VB edge lies ener¬ 
getically lower with respect to that obtained using the 
PBE functional. This causes several dye-related states to 
fall in the band gap just above VB. Therefore, the HSE 
functional gives a larger number of occupied gap states. 
In addition, the HSE functional corrects the band gap 
underestimation of the standard approximations to the 
exchange-correlation energy by admixing partial screened 
exact exchange. Therefore, unoccupied Ti 3d states are 
better described with the HSE functional leading to a 
significant shift of the CB up to higher energies. 

The cyanidin dye can form single and double bonds 
through its tail oxygens with the surface Ti ions on 
anatase. The latter is a chemical binding where the 
nearest surface oxygen captures an H from the tail OH 
group of cyanidin dye. Both adsorption modes bring sev¬ 
eral new occupied states above the VB and an empty 
state below the CB of Ti02 at both the PBE and the 
HSE levels of theory. The lowest lying transition start¬ 
ing from the HOMO-like dye-related gap state to the 
empty LUMO-like state below the CB takes place on 
the dye itself giving a sharp absorption peak in the vis¬ 
ible range.This can be understood as a factor limiting a 
subsequent charge injection into the nanowire and might 
end up with recombination of photoexcited electron-hole 
pairs. TA-St-CA shows mono- and bidentate binding 
modes on both of the anatase nanowires. When ad¬ 
sorbed, its intramolecular charge transfer character gets 
favorably modified toward a charge injection into the ox¬ 
ide. This might help reduce the recombination rates of 
charge carriers. The isolated hlled gap states originating 
from the bidentate binding of TA-St-CA sensitizer signif¬ 
icantly narrow the energy gap and effectively functional¬ 
ize the anatase nanowires to actively absorb a broader 
range of the visible spectrum. Inclusion of nonlinear sol¬ 
vation effects indicates dissociation of the monodentate 
TA-St-CA from the nanowires in water where the biden¬ 
tate TA-St-CA develops a moderate binding. In a more 
polar ionic solution, TA-St-CA needs a firmer anchor¬ 
ing to the oxide. Although the cyanidin molecule has 
a strong binding on the (101) and (001) nanowires even 
in an electrolyte, the absorption properties are weaker 
and might suffer from recombination of photogenerated 
electron-hole pairs. Bidentate TA-St-CA-|-nanowire sys¬ 
tems can achieve directional charge transfer excitation to 
increase charge injection probabilities, allow absorption 
in the wide range of visible spectrum with enhanced light 
harvesting, and exhibit moderate binding in solution re¬ 
quired to reduce degradation of possible device operation. 
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